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i— I Abstract 

o 

p<| We give an account, at nonexpert and quantitative level, of physics behind the CMB tem- 

perature anisotropy and polarization and their peculiar features. We discuss, in particular, how 
cosmological parameters are determined from the CMB measurements and their combinations with 
other observations. We emphasize that CMB is the major source of information on the primor- 

t-H dial density perturbations and, possibly, gravitational waves, and discuss the implication for our 

understanding of the extremely early Universe. 
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1 Introduction 

The existence of cosmic microwave background radiation (CMB) is one of the keystones of cosmol- 
ogy. It was discovered in 1964 by Pensias and Wilson, and it earned them the Nobel prize. CMB comes 
from the epoch when the Universe was much hotter and much younger than today. The approximate 
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age of the Universe at the time CMB was released equals 380 thousand years, while the current age 
of the Universe is about 14 billion years. Presently, the radiation has the temperature T = 2.72bK. 
Among all signals we receive directly, CMB radiation is believed to be one of the youngest)]] 

Later, the anisotropy of CMB was discovered. This anisotropy is very small. The largest is the 
dipole component of order 5T/T ~ 1CP 3 . It is believed to arise due to the motion of the Solar 
system with respect to CMB. All other multipoles are much lower, the maximal ones are of order 
5T/T ~ 10" 4 - 1(T 5 (see FigQ. 




-200 T(nK) • ?00 



Figure 1: The CMB temperature map, obtained by WMAP experiment [T]. The brighter a region is, 
the hotter radiation comes from it. 

It turns out that one can determine many parameters of our Universe by studying this anisotropy. 
CMB anisotropy is caused both by inhomogeneities in the early Universe and by the effects of propa- 
gation of CMB photons through the Universe at later time. An example of the latter is the Sunyaev- 
Zel'dovich effect, it is due to scattering of CMB photons off free electrons in interstellar plasma. This 
effect distorts CMB angular spectrum at small angular scales. There are other effects with similar con- 
sequences. For example, recombination is not an instantaneous process, the surface of last scattering 
has finite width, and this smooths out the CMB angular spectrum at small angular scales. Also, just 
before recombination the mean free path of photons drastically increases. The photons travel freely 
for rather long time, so they diffuse away from over-dense regions to regions of lower density, and this 
also smears out small inhomogeneities. This effect is called Silk damping. Because of these effects, 
the CMB anisotropy spectrum is not particularly informative from the cosmological point of view at 
angular scales of a few arc minutes and smaller. 

Last scattering of the CMB photons occurs almost at the same time as the recombination does, 
since the main process of interaction of photons with matter in this case is the Thomson scattering off 
free electrons. The temperature of last scattering depends on the cosmological parameters very weakly 
(weaker than logarithmically) , it is mainly determined by the potential of ionization of hydrogen and 
other world constants. The up-to-date values of the temperature and redshift of last scattering are: 

T r = 2970 K = 0.26 eV, (1) 
z r = 1090. (2) 

1 Earlier epoch of Big Bang Nucleosynthesis (starting from about 1 second after Big Bang) is probed by measuring 
light element abundances. 
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At the time of last scattering, the Universe was in matter-dominated regime, i. e., the contribution 
of mass density of nonrelativistic matter to the total energy density is the largest, or, in other words, 
the overall pressure is much smaller than the overall energy density: p tot "C ptot- 

CMB gives us the photographic picture of the Universe at the time of last scattering (modulo 
aforementioned effects). Hence, the CMB anisotropy reflects the fact that the Universe was slightly 
inhomogeneous at that time. The inhomogeneities in the following components of cosmic plasma 
contribute to the CMB anisotropy: 

1. Baryon-electron-photon plasma. 

2. Dark matter. 

3. Neutrinos. 

4. Gravitational waves. 

Rather complex structure of the CMB angular spectrum is due to the interplay between the effects 
caused by the first two of these components. 

The purpose of this review is to discuss what information about our Universe one extracts from 
the CMB data. There are two classes of properties one is interested in: 

1. Properties of primordial perturbations which, as we will see, were buil before the hot stage of the 
cosmological evolution. This is a key for revealing the origin of primordial perturbations and, 
ultimately, understanding the epoch preceding the hot stage. The existing data are consistent 
with the picture that primordial density perturbations and gravitational potentials induced by 
them had rather simple properties, consistent with the inflationary mechanism of their generation, 
and that there were no primordial gravitational waves. As we discuss in this review, future 
CMB measurements are sensitive to primordial gravity waves at the level predicted by simple 
inflationary models. The observation of the effects induced in CMB by primordial gravity waves 
would be an outstanding discovery. 

2. Properties of the Universe shortly before recombination and later, up until the present epoch. 
These include the energy balance (photons, neutrinos, baryons, dark matter, dark energy), spatial 
curvature, etc. 

Studying the CMB temperature anisotropy and polarization is sufficient by itself to measure some 
of these properties. There are numerous cases, however, when the CMB data alone are insufficient for 
determining some of the cosmological parameters, because these data are approximately degenerate in 
the parameters (i.e., they depend on certain combinations of the parameters). The degeneracy is lifted 
by invoking other cosmological data, as we illustrate in this review. 

There are good books and reviews on CMB, see, for example, |2l El SI El El El El 121 113 111] - The 
reader is invited to consult them, in particular, for references to original works. In this review we 
mostly refer to papers containing concrete results we make use of. 

2 Observables 

The temperature variation on celestial sphere, i. e., the function ST((f>,9) = ST(n), can be expanded 
in spherical harmonics: 



All observations today support the hypothesis that ai m are independent Gaussian random variables. 
Gaussianity means that 
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where P(ai m ) is the probability density for the random variable a\ m . For a hypothetical ensemble of 
Universes like ours, the mean values of products of the coefficients ai m would obey 
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We can measure only one Universe, but this formula is still used to extract the angular spectrum 
Ci from the data. For given I, there are (21 + 1) independent coefficients a/ m , so there exists an 
irreducible statistical uncertainty of order 5Ci jC[ ~ 1/ \/2l + 1, called cosmic variance. It is particularly 
pronounced at small I and, indeed, it is much larger than the experimental errors in this part of the 
angular spectrum (as an example, error bars in the left part of Fig. [2] are precisely due to the cosmic 
variance) . 

In the case of Gaussian random variables, all odd correlators are zero and all even correlators are 
expressed through the two-point correlator according to the Wick theorem. Under this distribution, 
the mean value of the temperature fluctuation squared equals: 
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Usually, graphs presenting the anisotropy spectrum show Di (for example, in Figj2|. It is worth noting 
that the approximate relationship between the multipole number I and angular scale is A9 ~ tt/1. 
Hence, the first peak in Fig jS] corresponds to temperature fluctuation of angular size of about 1°. 





6000 




5000 




4000 








3000 




2000 




1000 








■ WMAP 3 -year 
• ACBAR 

■ BOOMERANG03 




10 100 500 1000 1500 2000 2500 

t 



Figure 2: The angular spectrum of the CMB temperature anisotropy |12j . The line is a prediction of 
the standard KCDM model. The quantity in vertical axis is Di defined by @. 



Gaussianity of temperature fluctuations is a very interesting property. The linear evolution pre- 
serves Gaussianity, so this property means that the primordial density perturbations are Gaussian ran- 
dom field. We recall in this regard that vacuum fluctuations of free quantum fields are also Gaussian. 
Hence, we have an observational hint that the mechanism of the generation of primordial perturbations 
is the amplification of vacuum fluctuations of some linear quantum field(s). This is indeed the case in 
the inflationary scenario, and also in some other models. 
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3 Density perturbations and CMB temperature anisotropy 



3.1 Background metric 

To set the stage, let us begin with background. The background metric is the metric of homogeneous 
and isotropic space-time (Friedmann-Robertson- Walker metric): 



ds 2 = dt 2 - a{tf 



(/ ' J • r 2 {d6 2 ■ sin 5 
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where a(t) is the scale factor, k is the spatial curvature parameter, k = 0, +1, —1 correspond to the 
Euclidean 3-space, 3-sphere and 3-hyperboloid, respectively. In most part of this review we consider 
the case k = and use the equivalent form of the spatially flat metrics: 

ds 2 = dt 2 - a{t) 2 (dx\ + dx\ + dx\). (8) 

If we neglect the possible peculiar motion of cosmological observers, each observer has definite time- 
independent x-coordinates. Note that x\ + x\ + x\ is not a square of physical distance even in the 
spatially flat Universe (k = 0). The physical distance squared is a 2 (x 2 +X2+x 2 ), and this corresponds to 
the fact that the physical distance between two cosmological observers grows in time as a(t) increases. 
This is precisely the expansion of the Universe. The coordinates x are called comoving coordinates. 
In general, the dynamics of the scale factor is determined by the Friedmann equation, which reads: 



d\ 87rG k 
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(9) 



where p is the total energy density in the Universe. The energy balance equation (the relation dE = 
—pdV applied to the comoving volume a 3 ) is 

dp = -3(p + p)-, (10) 
a 

where p is pressure. For nonrelativistic matter (dark matter, baryons and also massive neutrinos at 
late epoch) one has p <^ p, and the energy density decreases as 1/a 3 (mass density gets diluted in 
inverse proportion to the comoving volume). For ultrarelativistic matter, for example, photons and 
neutrinos at early epoch, the equation of state is p = p/3 and p oc 1/a 4 (energy density decreases 
also because photons get redshifted). Equation ([9]) then implies that if the Universe is flat (k = 0) 
and filled with nonrelativistic matter, then the scale factor evolves as a oc i 2//3 (matter domination). 
If the Universe is filled with ultrarelativistic matter, then scale factor increases as a oc t 1 / 2 (radiation 
domination). There is also dark energy in our Universe, whose density does not change in time (or 
almost does not change in time), p\ = const. The Universe was radiation-dominated at early epoch, 
then there was long matter-dominated epoch, and recently dark energy has started to dominate the 
cosmological expansion. 

At the epochs we consider, the Friedmann equation is specified to 

TT n f a\ 2 8ttG , . k . 

H 2 = - = — Ow + PM + PA )--a= 11 
V a 1 6 a z 



8ttG 

Pc 



--' !( -) 4 + ^ M (?) 3 + fiA + ^(v) 2 



where the index "0" refers to the current epoch; /9 r ad>/°M)PA are energy densities of radiation, nonrela- 
tivistic matter (i. e., baryons and dark matter) and dark energy, respectively, p c = (8irG)/ (3Hq ), and 
^rad,M,A = /°rad,M,A/ ' Pc are the present energy densities in units of the critical density, while £l K is the 
relative present contribution of the spatial curvature to the current Universe expansion rate. Q's are 
relative contributions in the sense that f2 ra d + + + ^« = 1- 
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The coordinate t is the physical time, it is the proper time of cosmological observers. Sometimes 
another time coordinate is used, so-called conformal time 



In coordinates r\, x\, x%, X3 the metric ^ for k = reads: 

ds 2 = a{r 1 ) 2 (drj 2 -dx 2 ), (13) 

i. e., the metric is conformally flat. Since the light propagation satisfies ds 2 = 0, in coordinates 
77, xx, X2, X3 the light propagates as in the Minkowski space-time. 

3.2 The angular spectrum of the CMB temperature anisotropy 
3.2.1 Sachs Wolfe, Doppler and integrated Sachs Wolfe effects 

As we have already pointed out, the major source of the CMB anisotropy are the density perturbations 
and gravitational potentials associated with therrj^] In cosmological literature these are called scalar 
perturbations, since the relevant quantities are scalars under spatial rotations. In a certain gauge 
(conformal Newtonian gauge), metric with scalar perturbations is written as follows: 

ds 2 = (1 + 2<Z>)dt 2 - a 2 {t){\ + 2f)(ix 2 , (14) 

where $(x, t) and ^(x, t) are two gravitational potentials. In the Newtonian limit, coincides with 
the Newtonian potential. 

In the case of scalar perturbations, the observed temperature fluctuation of the CMB photons 
coming to us along the direction n from the point x where they have last scattered is given by the 
general formulcj^J 



6T 1 8p 1 



to 

(x) + $(x) + n- v 7 (x) + y dt ($-*), (15) 



where <5/? 7 (x) is the density perturbation of photons at last scattering, v 7 (x) is the velocity of the 
baryon-photon medium at that time, and the integral is evaluated along the photon world line from 
last scattering to detection. All terms here have transparent interpretation. The first term is the 
perturbation of the temperature of the photon gas at last scattering: in general p 7 oc T 7 and we have 
1 /4-(5 / o 7 (x) /p 7 = 5T 7 (x) /T 7 . The second term reflects the fact that photons get redshifted (blueshifted) 
when they escape from potential wells (humps). The first and second contributions together are called 
the Sachs- Wolfe effect. The third term is self-evident: this is the Doppler effect. Finally, the integral 



in (15) represents the integrated Sachs- Wolfe effect: as photons propagate through the Universe, they 
get red- or blueshifted due to time-dependent gravitational field. Namely, if a photon meets a potential 
well, it falls into it and gains energy. Suppose now that by the time the photon starts escaping from 
the well, the well has become shallower. Then the photon looses less energy when climbing the well, 
and the net effect is the blueshift of the photon. Since photons are relativistic, this effect is caused by 
both temporal (5goo oc <&) and spatial (Sgij oc components of metric perturbations. 

For scalar perturbations, both Doppler and integrated Sachs- Wolfe effects are numerically rather 
small. We simplify our discussion and consider the Sachs-Wolfe effect alone (the first two terms in 



(15). The main contributions here come from perturbations in the baryon-photon medium and dark 



2 We do not consider the effects due to structures in the relativeiy fate Universe, such as Sunyaev-Zel'dovich effect 
and gravitational lensing. Yet another possible source is gravity waves, see below. 

3 This formula does not account for aforementioned smearing effects, so it is valid, strictly speaking, only at relatively 
large angular scales corresponding to I < 1500. However, smearing leaves almost unaffected some properties of the CMB 
angular spectrum, such as positions of the peaks. We use this formula for qualitative analysis anyway. 
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matter. It is important that to a good approximation, baryons and photons make a single fluid before 
recombination; they are tightly coupled due to intense scattering of photons off free electrons and 
intense Coulomb interaction between electrons and protons^] 

Let us consider the effects due to baryon-photon medium and dark matter in turn. 



3.2.2 Oscillating part of the angular spectrum: baryon photon medium 

The oscillating part of the CMB temperature angular spectrum comes from sound waves in the baryon- 
photon medium before recombination. So, we have to understand the behavior of these waves in the 
expanding Universe. Instead of writing and then solving the equations for sound waves, let us consider 
a toy example, the case of inhomogeneities in massless scalar field. In fact, example is useful for 
demonstrating the main notions and ideas of treatment of various kinds of perturbations. The action 
for the scalar field is: 



S 



d^ 



The field equation thus reads: 



(fxdta 6 



d 
dt 



(a 3 0) + adidi<j) = 0, 
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1 
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(16) 
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l. e., 



< + m4) 



0. 



(18) 



Here, H denotes the Hubble parameter: H = a/ a. This equation is linear in <f> and homogeneous in 
space, so it is natural to represent <j) in terms of the Fourier expansion: 



x 



„ikx 



<p(k)d 3 k. 



(19) 



Clearly, the value of k for a given Fourier mode is a constant. However, k is not a physical wavenumber 
(physical momentum), since x is not a physical distance, k is called conformal momentum, while 
physical momentum equals q = 2n/\ = 2n/(a(t)Ax) = k/a(t). A here is the physical wavelength of 
perturbation, and it grows together with the expansion of the Universe. Accordingly, as the Universe 
expands, the physical momentum of a given mode decreases (gets redshifted), q(t) oc l/a(t). For a 



mode of given conformal momentum k, Eq. (18) gives: 



+ 3H0 + 



0. 



(20) 



This equation has two time-dependent parameters of the same dimension: k/a and H. One can 
consider two limiting cases: k/a <C H and k/a S> H. In the cosmological models with conventional 
equation of state of the dominant component (e. g., matter-dominated or radiation-dominated Uni- 
verse), one has H ~ 1/L, where L is the size of cosmological horizon. So, the regime k/a <C H is 
the regime when the physical wavelength is greater than the horizon size (super-horizon mode), while 
for k/a^> H the physical wavelength is much smaller than the horizon size (sub-horizon mode). The 
time when the wavelength of the mode coincides with the horizon size is called horizon crossing. This 
time is denoted by the symbol x. It is straightforward to see that at both radiation-dominated and 
matter-dominated epochs, the ratio k/(aH) grows. This means that every mode was at some early 
time super-horizon, and later it becomes sub-horizon. 

4 The tight coupling approximation is not valid for I > 1500 because of the Silk damping and other effects. Our 
discussion, however, remains valid at qualitative level even for large /. 
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For a super-horizon mode, we can neglect the term <j> ■ k 2 /a 2 in Eq.(18) in comparison with other 
terms. Then the field equation, e.g., in the radiation-dominated Universe (a oc t l / 2 ,H = l/2t), 
becomes 



The solution of this equation is 



+ |* = 0. (21) 



(t)=A + ^=. (22) 



This behavior is generic for all cosmological perturbations: there is a constant mode (A in our case) 
and a mode that decays in time. If we extrapolate the decaying mode Bj\ft back in time, we get very 
strong (infinite in the limit t — > 0) perturbation. For density perturbations this means that this mode 
corresponds to strongly inhomogeneous early Universe. Therefore, this mode is usually assumed to be 
absent for actual cosmological perturbations. Hence, for given k, the solution is determined by a single 
parameter, the initial amplitude A of the mode (p^. 

After entering the sub-horizon regime, the modes begin to oscillate. The equation in the sub-horizon 
regime is: 

0+^0+y<A = O, (23) 
and its general solution in the WKB-approximation reads: 

<K*) = -TTT cos ( / -7Ii\ dt> + ^0 I , (24) 




with the two constants being the amplitude A' and the phase ipQ . The amplitude A' of these oscillations 
is determined by the amplitude of the super-horizon initial perturbation, while the phase ipo of these 
oscillations is uniquely determined by the condition of the absence of the decaying mode, B = 0. 
Imposing these conditions yields 

0(i) = sin | / ~dt'\. (25) 




The decreasing amplitude of oscillations (p(t) oc l/a(t) and the particular phase xfjQ = —tt/2 in (24) 



are peculiar properties of the wave equation (18), as well as the radiation-dominated cosmological 



expansion. However, the fact that the phase of oscillations is uniquely determined by the requirement 
of the absence of the super-horizon decaying mode is generic. 

The perturbations in the baryon-photon medium before recombination — sound waves — behave 
in a rather similar way. Their evolution is as follows: 



p 1 



const, outside horizon, 

(26) 



const • cos yj v s - f J^Tjdt'^ , inside horizon, 



where v s = \J dp/ dp is the sound speed. The baryon-photon medium before recombination is almost 
relativisticj^] since p-Q < p 7 . Therefore, v s ~ l/v^- At the time of last scattering t r we have: 




5 This does not contradict the statement that the Universe is in matter-dominated regime at recombination. The 
dominating component at this stage is dark matter. 
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where 

t 



o 



is the comoving size of the sound horizon, while the physical size equals a(t r )r s . So, we see that 
the value of the density perturbation at recombination oscillates as a function of wavelength. The 
period of this oscillation is determined by r s . These oscillations are shown in Fig. [2] Indeed, one 
can show that gravitational potentials induced by the perturbations in the baryon-photon medium 
are negligible at recombination, while, as we mentioned above, the Doppler and integrated Sachs- 
Wolfe effects are also rather small. Hence, the main effect of the baryon-photon perturbations on 
CMB comes from the oscillatory first term in (15). The absolute values of perturbation^] in the 



baryon-photon plasma are maximal at conformal momenta k n obeying k n r s = im for integer n. 
Hence, their coordinate half wavelengths (distances between maxima of the absolute value) are equal 



to Ax n = ir/k n . Now, according to (13), photons propagate in conformal coordinates (ry,x) precisely 
like in the Minkowski space-time. Hence, an interval of coordinate length Ax at the time of last 
scattering is seen today at an angle Ax/rjQ, where rjo is the conformal time interval between the 
present and last scattering epochs, i. e., the comoving distance to the sphere of last scattering 
So, the strongest CMB temperature fluctuations have angular scales A9 n ~ r s /(nrjo). We recall the 
approximate correspondence A9 i — > it /I and conclude that the peaks in the CMB angular spectrum 
are expected to be at the positions 

l n ~ vrn— . (29) 
r s 

These peaks are indeed seen in Figj2] 

Let us pause here to record the approximate correspondence between the conformal momenta of 
perturbations and the multipole numbers of the CMB anisotropy they generate: 

I i — ► krjo- (30) 

This correspondence follows from the argument just presented and will be used later on. 

The very existence of the peaks in the CMB angular spectrum strongly suggests that the hot 
Big Bang stage of the cosmological expansion was preceded by another epoch with entirely different 
properties. That was the epoch when the cosmological perturbations were generated. Indeed, we 
have seen that the definite phases of sound waves which are behind the peaks are determined by the 
properties of the scalar perturbations in the early, super-horizon regime. So, the modes we observe were 
indeed super-horizon, their wavelengths exceeded considerably the size of the cosmological horizon. In 
the hot Big Bang theory, there is no causal connection between horizon-size regions. Hence, within this 
theory it is impossible to invent a mechanism which would generate these perturbations in a causal 
manner. We conclude that the generation of the cosmological perturbations must have occurred at 
some epoch preceding the hot Big Bang stage. Furthermore, the cosmological evolution at that epoch 
was such that the combination a(t)H(t) was small early on, and then it increased enormously towards 
the beginning of the hot Big Bang stage. Only in that case the physical momentum k/a(t) could 
exceed the Hubble parameter H at early times, the modes could be sub-horizon "in the beginning", 
and the perturbations could be generated by some causal mechanism. This is precisely what happens 
in the inflationary scenario. 

An alternative would be the generation of the density perturbations at the hot Big Bang stage in a 
causal way. This can only happen when a mode of perturbations is in the sub-horizon regime. In that 

6 According to ([H]), the relevant quantity is the temperature fluctuation squared. This is the reason for considering 
absolute values of density perturbations. 

7 This is true for intervals normal to the line of sight. One can show that the sound waves in the baryon-photon 
component, whose wave vectors are normal to the line of sight, indeed give dominant contribution to the oscillations 
in the CMB temperature angular spectrum. Contributions to the oscillating part of the angular spectrum coming from 
other waves get smeared out upon integration over the directions of the wave vectors. 
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Figure 3: The real CMB angular spectrum and inflationary prediction compared with the spectra 
obtained in models with casual mechanisms of the generation of density perturbations |13| . In all 
casually connected models (all models above except inflation), there are no oscillations. 



case, there is no reason for very definite phases of the sound waves in the baryon-photon component, 
and hence for oscillations in the CMB temperature angular spectrum. For example, Figj3] shows the 
angular spectrum in a model where density perturbations are generated by accretion on cosmic strings, 
and the spectra created by other casual mechanisms at the hot Big Bang epoch. All these spectra do 
not exhibit the oscillatory behavior. So, the generation of density perturbations at the hot Big Bang 
stage is ruled out by the CMB observations. 

The positions and strengths of the peaks are sensitive to various cosmological parameters, so from 
observations one can determine these parameters or at least their combinations. In particular, the peak 
positions depend on: 

1. The ratio of baryons to photons. 

2. The rate of the cosmological expansion, thus the density of dark energy and other forms of 
matter. 

3. The curvature of space. 



The size of the sound horizon entering (29) depends on the sound speed in the baryon-photon compo- 
nent. It is given by 

2 = dp = 1 d flI_Pl - 1 ^dT/T 1 1 

Vs dp 3^pb+^ 33^+4^ + { ' 

PB P-r Pi P"i 1 1 4 Pi 

The dependence of v$ and hence r s on /?b is rather weak, since Pb/p-/ "C 1 up to recombination. 
More sensitive to the baryon fraction are the heights of the peaks, see below. For given baryon 
abundance, precisely determined basically from the peak positions, the sound horizon serves as a 
standard ruler at recombination. The angle at which it is seen today depends on the expansion history 
after recombination and spatial curvature. 

The dependence on the expansion history of the Universe comes about in the following way. We 
know precisely the factor a(to)/a(t r ) determining how much the Universe has expanded from the time 
of last scattering (since a(to)/a(t r ) = T t /Tq and we know the present CMB temperature To and the 
temperature of last scattering T r ). The Universe could expand fast, and then the last scattering surface 
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would be close to us, or it could expand slowly, and the last scattering surface would be further away, 
so the angular scale at the peaks would be smaller. Notably, the expansion rate depends strongly on 
the value of the cosmological constant A (or, generally, on dark energy density); we discuss this point 
shortly in some detail. 

The dependence on curvature of space comes from the fact that the same interval viewed from the 
same distance is seen at different angles depending on curvature. For example, in the space of positive 
curvature (3-sphere) the angle is larger than the one in the Euclidean space, and the angle in the flat 
space is larger that the angle in the space of negative curvature (Lobachevsky geometry). So, both 
the curvature and the cosmological constant lead to similar shift of the peaks in the CMB angular 
spectrum, see Fig. [4] Therefore, there is a deg eneracy in the parameters of Q\ and see |14j . 
Actually, the degeneracy is lifted by making use of the results of cosmological observations unrelated 
to CMB. We give an example of the latter in Section [4] 




Figure 4: The effect of curvature (left panel) and of A (right panel) on the CMB temperature angular 
spectrum 



Let us illustrate the dependence of peak positions on the dark energy density by considering the 



spatially flat Universe. The Friedmann equation (12) in the flat Universe can be written as follows: 



8vrG 



PO (^rad(l + zf + fi M (l + Z? + «a) 



(32) 



where z = ao/a{t) — 1 is redshift. As we have already noted, the light propagates according to the 



equation ds = 0, i. e., dt 



a 2 dx z 



so 1] 



J (l/a)dt is the comoving distance traveled by light. To 



calculate the comoving distance which light has traveled since recombination, one uses Eq. (32) and 
writes 



Therefore, 



t 



dt 
a 



dz 



dz 



a 



da 



ao a 



dt 



a a 



-a H- 



jlt 



(33) 



dz 



a H{z) a HoJ y/n iad (l + zf + fi M (l + z) 3 + Q A 



(34) 



The effect on 770 due to is seen from this formula. Note that by definition f2 r ad + S1m + ^A = 1 in 
the spatially flat Universe, so the larger Q\ is, the larger is the integrand at z > 0. Since the positions 



of the peaks are approximately given by Eq. (29), these positions are shifted to the left at larger Oa- 
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3.2.3 Nonoscillating part of the spectrum: the effect of dark matter 

As we saw above, matter with finite sound speed produces oscillating angular spectrum of the CMB 
temperature anisotropy. But there exists another form of matter, which has zero (or nearly zero) sound 
speed. This is dark matter. Due to the fact that p = for dark matter, and hence v 2 = 0, there are no 
waves in this medium, and thus by itself it does not produce oscillations in the CMB angular spectrum. 
Nevertheless, dark matter density perturbations affect the CMB angular spectrum quite strongly, via 
the gravitational potentials they create. These gravitational potentials induce the CMB temperature 



anisotropy both directly, through the second term in (15), and through the induced perturbations in 
the baryon-photon medium, which falls into the potential wells. 

In the first place, we have to understand the properties of the gravitational potential <3?dm generated 
by dark matter perturbation. Let us consider a perturbation whose wavelength is much smaller than 
the horizon size at recombination. Then the Newtonian approximation is valid, and the gravitational 
potential is related to the mass density perturbation by the Poisson equation. We take into account the 
fact that the physical Laplacian written in terms of comoving coordinates is A = a~ 2 <9 2 /<9x 2 and write 
the Poisson equation in the expanding Universe in the following way (we sometimes omit subscript 
DM in $dm in this Section): 

k 2 

$ = 4ttG5 PbM - (35) 



a 2 (t) 

In matter-dominated regime <5/?dm//9dm oc a(t) for sub-horizon modes, while pr>M oc 1/a 3 . Therefore, 
5pDM oc 1 /a 2 and <£> is independent of time. Using the assumption of the flat spectrum of primordial 
perturbations, for which <5/»dm/pdm depends on k weakly (see below), we obtain: 

$DM(fc) OC p. (36) 

In fact, this relation is valid modulo logarithm, whose origin is the evolution of the dark matter 
perturbations at the radiation-dominated epoch; we are not going to elaborate on this point. 

To evaluate the effect of dark matter on CMB, we have to obtain the photon density perturbation 
Sp^/p^, which is induced by the gravitational potential ^dm- 111 this way we are going to find the 



first term in (15); as we have already pointed out, the third and fourth terms are rather small, and 
we neglect them for the sake of argument. To this end, let us consider the baryon-photon medium in 
static background gravitational potential. Baryon-photon plasma is then in the state of hydrostatic 
equilibrium. In Newtonian mechanics, the equation of hydrostatic equilibrium is: 

= -pV$. (37) 

In the case of relativistic medium, when p is of order p (but <I> is small, $ « 1), this equation is 
generalized to 

Vp = -(p + p)V$, (38) 

or, in the Fourier space: 

5p = -(p + p)$. (39) 

Here p and p are pressure and energy density of the baryon-photon medium before recombination. 
Using the equation of state p 7 = /3 7 /3, pb = we obtain: 

and, finally, 

ldp ^ = -(l + 3pB -)<S>. (41) 
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Let us pause here to point out that there is a subtlety in this result. Namely, suppose for the 
moment that baryon density is negligible, and consider static flat background space-time. Let there 
be dark matter which does not interact with photons, but creates the gravitational potential. Then in 
thermodynamic equilibrium the temperature of photons should be the same everywhere in space. This 
is in apparent contradiction with the result ( |41| , since the left hand side is precisely the perturbation 
of the photon temperature: bp^j p^ = 4<5T 7 /T 7 in view of the general property p 7 oc Tir. The reader is 



invited to reconcile thermal equilibrium with the result (41 ). 



temperature — the first two terms in (15) — is 



Coming back to CMB, we find from Eq. (41) that the Sachs-Wolfe effect on the observed CMB 



T 4p 7 

Notably, the net effect of dark matter is proportional to p-Q\ if there were no baryons, the gravitational 
potential would not affect CMB. 



According to (30) and (36), the dark matter contribution to the CMB angular spectrum is a smooth 
function of / that decays as 1/Z 2 . To large extent this contribution determines the overall decline of the 
spectrum, with oscillations superimposed on it, which is clearly seen in Fig. [2] (our Newtonian analysis 
is not valid for perturbations that remain super-horizon at recombination, so our discussion applies to 
the region I > 100 only). 

The effect of dark matter on CMB is so strong that the existence of dark matter could be deduced 
from CMB observations even in the absence of astrophysical evidence. As a digression, we point out 
that the presence of dark matter is crucial also for structure formation in the Universe. Indeed, at 
matter-dominated epoch, sub-horizon density perturbations in nonrelativistic matter (e. g., in baryons) 
grow as 5p/p oc a. At recombination, inhomogeneities in the baryon component were of order 10~ 5 , 
the Universe since then has the expansion factor a{to)/a{t r ) = z r + 1 ~ 1000, so the perturbations at 
present would be of order 10 -2 . They would not have left the linear regime of evolution, and thus would 
not form galaxies, stars, planets, men, etc. What saves the day is the property that the perturbations 
in dark matter grow well before recombination, and soon after recombination baryons catch up with 
them, as shown in Fig. [5] 

3.2.4 Interference 

The sum of the contributions of the baryon-photon and dark matter perturbations to the CMB tem- 
perature fluctuation has the following form: 

5 l = A-cos(kr s )-^$ DM . (43) 
T 4/9 7 

The first term comes from the baryon-photon acoustic oscillations, and the second one from dark 
matter perturbations. According to ([5]), the observables are quadratic in 5T/T, so an important issue 
is whether these two terms add up coherently. If they are coherent, we have: 

= A 2 cos 2 (kr s ) - l^cos(kr s )A$DM + ■ (44) 

TJ 2p 7 V 4 P 7 / 

Suppose now that the relative sign of A and $dm is positive. Then the peaks corresponding to 



cos kr s = — 1 (odd n in (29)) are enhanced, and the peaks corresponding to cos kr s = 1 (even n) are 
suppressed due to the second, interference term. Since $ oc 1/A; 2 , the effect is more pronounced for 
the first peaks. And this is exactly what is observed. 

Let us recall that the density perturbations are random fields. Hence, the interference is possible 



only if the baryon-photon perturbations (the coefficient A in (43)) and dark matter perturbations 



($DM in (43)) are proportional to one and the same random field, call it 1Z(\t). This means that the 
baryon-photon and dark matter perturbations are of one and the same origin, with 1Z being the common 
amplitude of primordial scalar perturbation. In fact, CMB data are consistent with the property that 
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Radiation domination 



Matte 



r domination 



A domination 




Figure 5: Schematic plot of the evolution of density pertirbations in different components. Here, 
5 = Sp/p, and $ is the gravitational potential. The left dashed vertical line is the time of horizon 
crossing of a mode considered. 



primordial scalar perturbations are in the adiabatic mode. The definition of the adiabatic mode is 
that the particle content is one and the same throughout the Universe. In other words, adiabatic 
mode would appear if one contracts or expands some regions of the Universe without changing the 
chemical composition of matter in these regions. The invariant and time independent characteristic of 
the baryon abundance is the ratio ub/s of the baryon number density to the entropy density. Hence, 
in the adiabatic mode hb/s is constant in time and space. Likewise, the ratio of the number density of 
dark matter particles to the entropy density udm/s is a universal constant. Since s oc T 3 and 
for the adiabatic mode in super-horizon regime we havd^j 



8ppM _ 3 Spy _ 6pB _ 3 Spu 
Pdm 4 p 7 p B 4 p v 



K(k). (45) 



After the horizon entry, the perturbations in the baryon-photon and dark matter components start 
to evolve. This evolution is linear (perturbations are small) up until recombination, so we have at 
recombination 

$ DM (k) = T DM (k)K(k), (46) 
A(k) = r B7 (fc)ft(k). (47) 

The functions Tdm and Tb 7 are called transfer functions. They describe how the perturbations in 
different media evolve. Yet $dm an d A at recombination are proportional to one and the same random 
field 1Z, and they add up coherently. 



The relations (451 are valid in the conformal Newtonian gauge. A convenient gauge-invariant definition of TZ is the 



spatial curvature of comoving hypersurfaces. 
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In principle, besides the adiabatic mode, there could exist so-called CDM entropy mode and baryon 
entropy mode. In the CDM entropy mode, there are no primordial perturbations in the baryon-photon 
medium, perturbations exist only in dark matter. Such a situation persists until the horizon crossing, 
due to causality. The inhomogeneities in the baryon-photon medium are generated at the time of 
horizon crossing, and afterwards they oscillate as follows: 



oc sin 



V 



/' 



k 



•a(t>] 



dt' 



(48) 



We emphasize that the phase of oscillations is different for adiabatic perturbations (cosine rather than 
sine, see (27)). So, if primordial perturbations were in the CDM entropy mode, the peaks in the 
CMB angular spectrum would be shifted, see Fig. [6] Exactly the same result holds for the baryon 
entropy mode. Presently, the admixture of the CDM entropy perturbations is ruled out by the CMB 
observations at the level of about 10%. 
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Figure 6: The effects of adiabatic and entropy perturbations on the CMB angular spectrum [IT] . 



Purely adiabatic mode of perturbations is very natural. If our Universe was in complete thermal 
equilibrium (including chemical equilibrium) at some early stage of the hot Big Bang epoch, the 
only mode of perturbations is adiabatic. Physical processes that generated dark matter and baryon 
asymmetry, if that happened at the hot Big Bang epoch, were the same everywhere in space, so the 
values udm/s and hb/s are the same everywhere^] Reversing the argument, if there exists an entropy 
mode in our Universe, the generation of the baryon asymmetry and/or dark matter must have happened 
before the hot Big Bang epoch. Hence, the detection of any of the entropy modes would be a strong 
signal for rather unconventional cosmology. 

To end up with the interference, we note that the effect is proportional to p-Q and hence Ob- The 



same holds for the third, nonoscillating term in (44), which gives particularly strong contribution in 
the region of the first peak. 

In this way the CMB temperature anisotropy is very sensitive to ps/p-y and, therefore, to the 
baryon asymmetry parameter 7/b = n^/n^ and Qr. This is shown in Fig|7] The CMB result for the 
latter is: 

n B = 0.045 ± 0.003. (49) 

Notably, it is consistent with the determination of t?b from Big Bang Nucleosynthesis. 

9 The argument is, in fact, even more general: the adiabatic mode is the only mode also in the Universe where particles 
of all types are decay products of one and the same field, like inflaton in the inflationary scenario. 
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4 Baryon acoustic oscillations 

We now digress from the discussion of CMB and introduce one of the independent methods of measuring 
combinations of cosmological parameters. This method has to do with a peculiar property of matter 
density distribution in the Universe. 

At the time right after recombination, the total matter density contrast (dark matter and baryons) 

is 

— = 1 = n + 7S ' ( 50 ) 

P Pb + Pdm »m Pb s 'm Pdm 

As we have seen above, right after recombination 5ps oscillates as a function of k, while Sp-^M is a 
smooth function of k. After recombination, the matter density distribution is essentially constant in 
time (in comoving coordinates), because baryons have decoupled from photons, and the sound speed 
in the baryon component is essentially zero. Therefore, the power spectrum of the total matter density 
perturbations today has an oscillating in k part. This phenomenon is called baryon acoustic oscillations 



(BAO). In the adiabatic mode, the two terms in (50) add up coherently, so there is an interference 
between them in (5p/p) 2 . The interference is important, since if it did not exist, the oscillations would 
be the second order effect in the small parameter S7b /^m- D ue to the interference, the oscillations are 
of the first order in J)b/^m : 

*P) a = ^B(k) + ^C(k)cos(kr s ). (51) 

P J I'M "M 

This effect has been detected in the analysis of the distribution of galaxies in the Universe at rather 
small redshifts, z ~ 0.2 and z ~ 0.35, see |15j . 

There is a simple interpretation of the effect. In the adiabatic mode, the overdensities in the 
baryon-photon medium and in the dark matter are at the same place before horizon crossing. But 
before recombination the sound speed in baryon-photon plasma is of the order of the speed of light, 
while the sound speed in dark matter is basically zero. So, the overdensity in baryons generates the 
outgoing density wave after horizon crossing. This wave propagates until recombination, and then 
freezes out. On the other hand, the overdensity in the dark matter remains in its original place. The 
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Figure 8: The correlation function of matter density versus distance [16J. (1) - The correlation 
function in the model with baryons and CDM, (2) - the correlation function in pure CD M-model 
without baryons. 



current distance from the overdensity in dark matter to the front of the baryon density wave equals- °| 
150 Mpc. Hence, there is an enhanced correlation between matter perturbations at this distance scale, 
which shows up as a feature in the correlation function, see Fig. [8] In the Fourier space, this feature 
produces oscillations |15) . 

Since the sound speed in the baryon-photon medium before recombination is well known, these 
150 Mpc is a robust prediction, so BAO provides a standard ruler at low redshifts. In principle, it can 
be used to determine the Hubble parameter as function of z at low z (if the ruler is directed along the 
line of sight, its absolute length I is related to its measured extension in redshift by Az = H(z)l) and 
also its angular diameter distance D a (z) (if the ruler is normal to the line of sight, the angle at which it 
is seen is A9 = 1/ D a (z)). The latter again depends on the expansion history (now at low z) and spatial 
curvature. In practice, one determines from observations a combination [D^(z) / 'H(z)] 1 / 3 , which is a 
rather complicated function of Q\ and f2 K . Importantly, this is a different combination as compared 
to the one best measured in the CMB observations. Hence, BAO are used to lift the degeneracy in 
parameters inherent in the CMB data, see |14| . 



5 Primordial perturbations 

5.1 Scalar power spectrum and scalar tilt 



We have seen in Section |3.2.2 that the oscillations in the CMB angular spectrum can be naturally 



explained if we assume that the scalar perturbations existed already at the beginning of the hot Big 
Bang epoch of the cosmological expansion. Also we mentioned that the scalar perturbations are 
observed to be Gaussian, and thus all information about them is contained in the two-point correlation 

10 We remind that the separation at Figgis given in h' 1 Mpc, where h = if/100 km/(s ■ Mpc), H is a current Hubble 
parameter. Currently h « 0.7, so 100 ft -1 Mpc roughly corresponds to 150 Mpc. 
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function: 



(K(k)K^k')) = ^5(k-k>), (52) 

where lZ(k) parametrizes the value of the initial perturbation, and P{k) is called power spectrum. In 
the isotropic situation, P(k) depends only on the magnitude of vector k, but not on its direction. The 
mean square of fluctuation is 

oo 

(7£ 2 (x)) = ( j e lkx R(k)d 3 k J e- ik ' x R*{k')d 3 k'} = J d 3 k^^ = J ^{k), (53) 

o 

where V{k) = k 3 P(k)/(2ir 2 ) is also called power spectrum. The quantity V{k) measures the contribu- 
tion of a logarithmic interval of k to the total power of fluctuations. Usually, V(k) is approximated as 
follows, 



k x ns ~ X 



V s {k) = A s \ — \ , (54) 

where k* is some fiducial momentum, A s is the amplitude of the power spectrum at this momentum, 
and n s is called the spectral tilt. Subscript s refers to scalar perturbations. In 1960's, Zel'dovich and 
Harrison predicted the flat spectrum of perturbations, that is the spectrum with n s = 1. Approximately 
flat, Harrison— Zel'dovich spectrum is consistent with the inflationary scenario. 

CMB data (and also data on the distribution of galaxies in space) are obviously sensitive to both the 
amplitude A s and the tilt n s . The amplitude determines the overall magnitude of the CMB temperature 
anisotropy, while the tilt gives the relative power at high and low multipoles. Both parameters can be 
measured rather precisely. The current values with k* = (500 Mpc) -1 are |14| : 

A s = (2.46 ±0.09) • 10" 9 , (55) 
n s = 0.960 ±0.014. (56) 

Note that the Harrison-Zel'dovich value n s = 1 is somewhat disfavored. The quoted values are 
obtained, however, under the assumption of the absence of tensor perturbations (see below) and exact 



form (54) of the primordial power spectrum. Relaxing these assumptions makes the allowed intervals 
of A s and n s wider. Note also that A s characterizes the perturbation squared; the amplitude of 
perturbation itself is of order 

V^-5-10- 5 . (57) 



According to the formula (45 ), IZ(k) is a relative amplitude of initial perturbations. Also, \fA~ s measures 



H(k*), so \/Ag is approximately equal to the relative CMB temperature fluctuation 5T/T. 

5.2 Inflation and generation of perturbations 

Inflation is fast, nearly exponential expansion of the Universe: 



a(t) oc exp Hdt^j , H m const. 



(58) 



Inflation occurs if the Universe is filled with a scalar field ip (inflaton) which has nonvanishing scalar 
potential V (up). The homogeneous field <p satisfies the equation 

<p + M<p = - d --. (59) 
dip 
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If H is large (the Universe expands rapidly) and dV/dip is small, the field <p varies slowly in time. The 
Friedmann equation in this case is H 2 = 8irGV(tp)/3, so if (p varies slowly, then V(<p) and thus H also 
vary slowly. This gives a self-consistent mechanism of inflation. 

The perturbations about the homogeneous field ip evolve according to the same equation as in 



Section 3.2.2| (one can show that the effect of V(<p) on 5(p is negligible): 



k 2 

5<p + SHdip + -zSip = 0. (60) 
a z 

Recall that if k/(aH) 3> 1, the mode is sub-horizon, and if k/(aH) <C 1 the mode is super-horizon. In 
the inflationary regime one has a(t) oc exp(Ht), H is almost constant, and k/(aH) decreases in time. A 
mode of a fixed k is sub-horizon at the beginning of inflation. Then it exits the horizon, and the mode's 
evolution terminates. Note that this sequence of events is reverse as compared to what happens in 
the matter-dominated or radiation-dominated Universe. So, to find the perturbations in the Universe 
towards the end of inflation, one needs to calculate the inflaton perturbation at the moment of the 
horizon crossing^] 5tp x ■ There are two quantities of the right dimension: k/a and H. But at the time 
of the horizon exit one has k/a = H. So, by dimensional argument, one obtains that at the time of 
the horizon exit 

6<p* ~ H, (61) 

and 5ip x is almost independent of k. This value remains intact until the end of inflation. 

After inflation ends, the energy stored in the inflaton field gets converted into heat, and the Universe 
enters the hot Big Bang epoch. The inflaton perturbations get reprocessed into scalar perturbations 
(energy density perturbations of the hot gas of particles and gravitational potentials). These are 
automatically in the adiabatic mode, since the perturbations in all forms of matter have one and the 



same origin everywhere in space 12 Since the parameters in the inflating Universe are almost time- 
independent in the time interval when all interesting modes exit the horizon, inflation produces nearly 
Harrison-Zel'dovich, flat spectrum (we have already noticed that 5(p x is almost independent of k). Yet 
the parameters are not exactly time-independent at inflation, so the predicted spectral tilt (n s — 1) is 
small but nonzero. It can be positive or negative, depending on the shape of the scalar potential V(<p). 
In particular, it is negative for the simplest power-law potentials like 

U(^) = ^V or V(p) = -^ A . (62) 
Finally, the primordial perturbations in the field cp have to be small in amplitude to generate small 



amplitude of the scalar perturbations (57). Thus, 5ip can be treated in the linear approximation. Linear 
quantum fluctuations of scalar fields are known to be Gaussian random variables, so inflation naturally 
explains this property of density perturbations. We conclude that the inflationary mechanism of the 
generation of scalar perturbations is consistent with everything we know about them. 

Actually, the derivation of approximately flat spectrum does not depend on whether we deal with 
scalar or tensor field. So, inflation generates also tensor perturbations (transverse traceless perturba- 
tions of spatial metric htj, i. e., gravitational waves). We obtain the same picture for them: primordial 
tensor perturbations are Gaussian random field with almost flat power spectrum. The standard approx- 
imation for primordial tensor perturbations spectrum (omitting the details concerning polarization) 
is: 

V T (k) = A T f . (63) 



3.2.2 



A 



11 Note that the meaning of the moments denoted by the symbol x is different in this Section and in Section 
given mode crosses the horizon at least twice, for the first time at inflation (it turns from sub-horizon to super- horizon) , 
and for the second time after inflation (it becomes sub-horizon again). In this Section x denotes the first crossing, and 
in Section 



3.2.2 



the second one. 



12 In inflationary models with more than one scalar field, the generation of an admixture of entropy modes is also 
possible. 
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Note the flat spectrum corresponds to tlt = 0, unlike in (54). It is convenient to introduce the 
parameter r = Vt/Vs which measures the ratio of tensor to scalar perturbations. The simplest theories 
of inflation like (62 ) predict r ~ 0.1 — 0.3. The effect of tensor perturbations has not yet been reliably 



observed, but there is great effort to discover it. The allowed regions in the plane (n s ,r) are shown in 
Figj9] The detection of tensor perturbations is an extremely important test for the simplest inflationary 
models. The flat primordial spectrum of scalar perturbations is predicted by some theories other than 
inflation, but at the moment only inflation predicts the flat spectrum of primordial gravitational waves. 
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Figure 9: The allowed regions of n s and r [14] . The contours show 68% and 95% CL from 
WMAP+BAO+SN. The label "HZ" corresponds to the flat, Harrison-Zel'dovich scalar spectrum in 
the absence of tensor perturbations. Also shown are the predictions of various inflationary models, 
whose potentials are given on the right. iV-flation is a model with numerous inflaton fields. 



The tensor perturbations behave at the hot Big Bang epoch exactly in the same way as massless 



scalar field studied in Section 3.2.2 These are purely metric perturbations (essentially no perturbations 
in matter): 

ds 2 = dt 2 - a 2 (t)[Sij - hij{x, t)]dx l dx j , (64) 

where hi* is a transverse traceless tensor, dihy = 0; ha = 0. Before the horizon re-entry, metric 
perturbation hij stays constant in time; after the horizon entry it oscillates and its amplitude decreases 



as \ ja. For modes that cross the horizon at radiation domination one has (compare with (25 ); we omit 
indices 



^l-sinU ! 
a(t) Sm \ k J a(i') 



h(k,t) = h (l) (k) ■ ^ -sin [ k I ^ | , (65) 



where hn\ is the primordial amplitude (Gaussian random field). The phase of oscillations is different 
for modes that cross the horizon at matter domination. 

For nearly flat primordial power spectrum, the shorter the mode, the earlier it crosses the horizon, 
the smaller a(t x ), the lower the present amplitude. For this reason it is almost hopeless to detect relic 
gravity waves directly by existing and planned detectors of gravitational waves. Presently, the best 
constraints on tensor perturbations come from the analysis of the CMB temperature anisotropy. The 
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effect of gravitational waves on the CMB temperature anisotropy is given by: 

to 

dt hijrfn 3 , 

i, 



6T 



(66) 



where n l is the unit vector along the photon path and the integration is performed along the photon 
world line. This is nothing but the integrated Sachs- Wolfe effect; other effects we considered in 
Section [3. 2.1 are absent for tensor perturbations. The modes which are super-horizon at recombination 
(I < 50) have not changed since inflation, as they have never been sub-horizon after inflation but before 
recombination. So, their spectrum at recombination is flat, and they generate flat contribution to the 
CMB temperature angular spectrum at low I. At higher multipoles, shorter gravity waves contribute 



recall the relationship (30)), their amplitudes decrease with k, and hence their effect on the CMB 



temperature decreases with This is shown in Fig. 10 The oscillations in the angular spectrum are 



due to the oscillatory behavior in ( 65 ) . 



Clearly, it is hard to disentangle the contributions of scalar and tensor perturbations to the CMB 
temperature anisotropy. In particular, tensor perturbations and red scalar tilt (n s < 1) have similar 
effect of the enhancement of low multipoles (this is the basic reason for the shape of the allowed region 
in Fig. [9]). Potentially more powerful way to detect primordial gravity waves is the measurement of 
the CMB polarization. 




Figure 10: The effect of gravity waves on CMB [llj. 



6 CMB polarization and primordial gravitational waves 
6.1 The physical mechanism of polarization 

CMB polarization is generated by a simple mechanism: light gets polarized when it scatters off free 
electrons. Indeed, when the light is scattered off an electric charge, the electromagnetic wave induces the 
oscillations of the charge in the direction of electric field. Electric field is perpendicular to the direction 
of the wave propagation. The oscillating charge emits radiation only in directions perpendicular to 
the direction of its oscillations. So, if the initial light is polarized in the plane formed by the direction 
of initial propagation and direction of scattered light (scattering plane), smaller amount of radiation 
is emitted as compared with incoming light polarized normal to the scattering plane. A photon with 
polarization normal to the scattering plane scatters at the cross-section 

da 3(Jt 



dfl 8vr 



(67) 
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and its polarization remains unchanged. If the initial polarization is in the scattering plane, it remains 
in that plane, and the cross-section is 



da 3a T 2q ( RQ \ 
—— = cos 0, 68) 

where 9 is the scattering angle. If at some position the flux of photons is anisotropic, the radiation 
scattered there towards an observer is linearly polarized. The local anisotropy before the last scattering 
can be caused by both temperature inhomogeneities and gravitational waves. The effect involves the 
factor A 7 /A(/c), where A 7 is the photon mean free path just before last scattering, and A(A;) is the 
physical wavelength of the perturbation which is responsible for the local anisotropy. So, the estimate 
for the degree of polarization P induced by temperature inhomogeneities is 

A(jfe) T ' [W) 

For I ~ 100 (this roughly corresponds to the horizon at recombination) the degree of polarization is 
approximately 10 -6 . 

6.2 E— and B— modes 

Consider a polarimeter oriented along a vector n normal to the line of sight. The intensity it measures 
is given by 

I n = ((En) 2 ) = n a n b (E a E b ) = (E 2 )P ab n a n b , (70) 
where the polarization tensor is defined as follows: 

_ (E a E b ) 

p ab-^^y (71) 

All averages here are the time averages over the period of the wave. For linear polarization, the tensor 
P ab is real. Additionally, it is symmetric and satisfies the constraint Tr(P) = P aa = 1. A symmetric 2- 
tensor has three parameters, and since the polarization tensor satisfies one additional trace constraint, 
it is characterized by only two parameters. The eigenvectors of this tensor are the directions of maximal 
and minimal intensity. As is clear from the above discussion, the direction of maximal intensity near a 
hot spot is perpendicular to the direction towards the spot. Near a cold spot the direction of maximal 
intensity is the direction towards the spot. 

CMB polarization can be visualized by showing the direction of maximal polarization and the 
degree of polarization V = \/\ — 4detP. The actual CMB polarization map is shown in Fig 11 

It is convenient to parametrize the polarization tensor as follows: 

Pab ~ \5ab = -{VaVb}P E ~ KV 6 V C }P B , (72) 

where {^4 a b} = {A ab + A ba — 5 ab A cc )/2. The modes generated by the potentials Pg and Pg are called 
P-mode and P-mode, respectively, in distant analogy with electric and magnetic fields. The P-mode is 
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obviously parity-odd, while the P-mode is parity-even. Scalar perturbations generate only E-modt 
while tensor perturbations generate both modes. This is because that parity-odd B-mode can noy 
generate parity-even scalar perturbations. Similarly to the expansion of temperature, one expands 
potentials of CMB polarization in spherical harmonics; omitting Z-dependent factors, we write 

P EiB (n) = J2^t BY Un). (73) 

Lm 



13 Some admixture of B-mode is generated at small angular scales from _E-mode via gravitational lensing by structures 
(clusters of galaxies, etc.) in the late Universe. 
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Figure 11: CMB polarization map in WMAP K band. [I] 



Equivalently, a lr ^ are the expansion coefficients of P a b in parity-even and parity-odd tensor spherical 
harmonics. Polarization observables are the correlators 

C EE (l) oc (af m afj, (74) 
C BB (l) oc (a B m a B J. (75) 

Another possible correlator, (aj^aj^), vanishes because of its odd parity (a E are parity-even and a B 
are parity-odd). Since CMB polarization is caused predominantly by local anisotropics in the medium, 
which are generated by the temperature inhomogeneities, one naturally expects that polarization will 
correlate with the CMB temperature fluctuations. Hence, one introduces the cross-correlator 

C te (1)k (a E m afj. (76) 

The cross-correlator C TB {1) oc (aj^af^ also vanishes because of the parity argument. In observations 
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only C TE and C EE correlators were discovered so far, and there exist upper limits on C BB , see Fig 

Scalar perturbations contribute to polarization mainly because of the motion of the baryon-photon 
medium. Photons coming along the direction of motion are hotter than photons coming from the 
opposite direction. This is precisely the Doppler effect that produces local anisotropy before the very 
last scattering, and hence CMB polarization. It is clear from the discussion in Section |6.1| that this 
local anisotropy by itself is insufficient for producing CMB polarization (effects due to photons coming 
from opposite directions would cancel out), so there is still a suppression factor A 7 /A(fc). For the same 
reason, the effect of density inhomogeneities on CMB polarization is of the second order in A 7 /A(fc). 
In acoustic waves, the phase of oscillations of velocity is shifted by n/2 from the phase of oscillations 
of density. The effect of this shift is that the oscillations of Cj E are shifted by half-period with respect 
to oscillations in Cf T (recall that the oscillating part of the CMB temperature anisotropy spectrum 
Cf T comes mainly from 8p~, recall also the correspondence (30); in the previous Sections we did not 



write the superscript TT in C{): zeroes of Cf E should coincide with maxima and minima of Cf T 



This can be seen in Fig. 13 



The effect of gravitational waves on CMB is much more interesting. At relatively low multipoles, 
primordial gravitational waves (tensor perturbations) are the only possible source of considerable B- 
mode of CMB polarization. Thus, the detection of the B-mode will be the detection of primordial 
gravity waves. Measuring their power spectrum will be of paramount importance for understanding 
the origin of cosmo logical perturbations. 
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Figure 12: Correlators of various polarization potentials and temperature as functions of I |17| . C EE 
and Cf B are multiplied by Tq , and Cf is multiplied by To, so all these quantities have the same 
dimension as C TT . The BB-mode has not yet been detected; the lower right panel shows the existing 
limits and predicted spectrum generated via gravitational lensing by structures in the late Universe in 
a model without primordial gravitational waves. 



As we have already mentioned, the smaller the wavelength of a gravitational wave, the earlier it 
becomes sub-horizon and the smaller its amplitude at recombination. So, the effect of gravitational 
waves is small for small wavelengths, and hence for large I. On the other hand, the effect of gravitational 
waves at recombination contains the factor (A 7 /A(A;)) 2 , so it is small for very large wavelengths and 
thus for very small I. The maximal effect is at I ~ 50. We note, however, that the interaction of 
gravitational waves with light, which leads to CMB polarization, occurs not only at the time of last 
scattering. Approximately at z ~ 10, intergalactic gas was reionized by the light of first very massive 
stars with M ~ 100 M®. Tensor modes of wavelengths approximately equal to the size of the horizon 
at this reionization epoch are most important in this situation. This means that the effect is most 
pronounced at I < 10. This is seen in Fig. 14 , a wide peak at / ~ 2 — 7 is precisely due to reionization. 
The optical depth at reionization (i. e., the probability for a photon to be scattered) is about 0.1, this 
is the reason for the smallness of the effect. 

Special experiments intended to measure the CMB polarization, such as CMBPole, are being pre- 
pared. The sensitivity of CMBPole to r = At/A s is approximately 10%, and it is planned to increase 
the sensitivity to the level of 1% in the next generation experiments. 
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Figure 13: Effects of scalar perturbations: CT E compared with CT T [T|. Zeroes of Cf E approximately 
coincide with maxima and minima of CT T and vice versa. 




Figure 14: The effect of scalar perturbations (left panel) and gravity waves (right panel) on CMB [18J. 
Primordial scalar and tensor spectra are assumed to be flat. The relative normalization is arbitrary. 

7 Summary 

To summarize, measurements of the CMB temperature anisotropy and polarization is a powerful tool 
in cosmology. Combined with other methods, they are capable of determining with high precision the 
cosmological parameters characterizing the recent Universe and revealing the properties of the primor- 
dial perturbations. The current picture of the Universe is quite simple. As to the recent Universe, the 
data are consistent with the spatially flat ACDM model whose ingredients are time-independent dark 
energy density, cold dark matter, baryons, electrons, photons and fairly light neutrinos [m u < 0.2 eV 
for each of the neutrino species). Known properties of scalar perturbations are also simple: they 
were generated before the hot Big Bang epoch, have no decaying super-horizon modes, are adiabatic, 
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Gaussian and have flat or nearly flat power spectrum. Tensor perturbations have not been discovered, 
though the current bound on their amplitude is not particularly strong. All these properties are con- 
sistent with the inflationary mechanism of the generation of cosmological perturbations, but inflation 
is not the only option for the moment. 

What can one expect to be discovered in future? On the recent Universe side, the major issues 
are whether the dark energy density is exactly constant in time or not, and what are dark matter 
particles. CMB alone will not shed much light on these issues, but its analysis will be instrumental in 
combination with other methods. Concerning the cosmological perturbations, one intrigue is whether 
the scalar power spectrum is exactly flat or not. Even more interesting are tensor perturbations whose 
primordial amplitude is predicted by the simplest inflationary models to be quite large; furthermore, 
the smoking-gun prediction of inflation is nearly flat tensor spectrum. Unexpected discoveries cannot 
be excluded either, like sizeable non-Gaussianity or admixture of entropy modes in scalar perturbations. 
All these are the tasks for future CMB experiments, which will thus serve as windows to the extremely 
early cosmological epoch preceding the hot Big Bang stage. 
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